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The self-assembly of iron dots on the insulating surface of NaCl(OOl) is investigated experimen- 
tally and theoretically. Under proper growth conditions, nanometer-scale magnetic iron dots with 
remarkably narrow size distributions can be achieved in the absence of a wetting layer. Further- 
more, both the vertical and lateral sizes of the dots can be tuned with the iron dosage without 
introducing apparent size broadening, even though the clustering is clearly in the strong coarsening 
regime. These observations are interpreted using a phenomenological mean-field theory, in which a 
coverage-dependent optimal dot size is selected by strain-mediated dot-dot interactions. 



PACS numbers: 81.07.Ta,61.46.+w,68.37.Ps,68.35.Md 

Clustering on surfaces by nucleation and growth dur- 
ing atom deposition has been an important subject in ba- 
sic and applied science for decades gj. Recent efforts 
have been focused on searching for methods to obtain 
nanometer-scale clusters with narrow size distributions. 
Such clusters or quantum dots are potentially valuable 
for optical, electronic, and magnetic device applications, 
but mass production of such structures by lithography or 
etching-based fabrication has proved to be exceptionally 
challenging H |) . Alternatively, it has been realized that 
the strain energy associated with the lattice mismatch 
between the dot and the substrate materials can be ex- 
ploited to induce self-assembled formation of quantum 
dots with narrow size distributions. This has generated 
much excitement, particularly in the area of semicon- 
ductor quantum dots |j[ ||, ||, ||, 0, ||. In such cases, 
the growth of the dots often proceeds in the Stranski- 
Krastanow (SK) mode, which is characterized by the 
presence of a wetting layer prior to three-dimensional 
(3D) clustering. To date, the precise mechanism for 
size selection in semiconductor quantum dot systems re- 
mains a subject of active debate f§ ||, |[ [§: some at- 
tribute them to strain-induced thermodynamic equilib- 
rium states, while others associate them with metastable 
configurations due to kinetic limitations. 

Although improved size uniformity can be achieved in 
quantum dot formation via the SK growth mode, the 
presence of a wetting layer is often undesirable, particu- 
larly for electronic and magnetic device applications of 
metallic/magnetic quantum dots. For this reason, it 
is preferred to fabricate quantum dots in the Volmcr- 
Weber (VW) growth mode, which is characterized by 
immediate 3D clustering on the substrate surface. In- 
deed, considerable recent efforts have been devoted to 
metallic/magnetic quantum dot formation on various 



substrates |T^, [Tl], [T^] , but no significant size uniformity 
has been achieved in such studies. 

In this Letter we investigate the self-assembly of iron 
dots on NaCl(OOl), an insulating substrate, by thermal 
deposition and variable-temperature atomic force mi- 
croscopy in ultrahigh vacuum. We show that, by properly 
choosing the growth conditions, nanometer-scale mag- 
netic iron dots with remarkably narrow size distribu- 
tions can be achieved in the absence of a wetting layer 
(VW growth). Moreover, by changing the dosage of iron, 
we can tune both the vertical and lateral sizes of the 
dots without introducing apparent size broadening, even 
though the clustering is already in the strong coarsen- 
ing regime, signified by the decrease in dot density as a 
function of the iron dosage. The preserved narrowness in 
the island size distributions is in clear contradiction with 
the expectations of existing understanding of clustering 
on surfaces in the coarsening regime [^. We interpret 
these observations within a phenomenological mean-field 
theory, in which a coverage-dependent optimal dot size 
is selected by the competition between the self-energy 
of a dot, i.e. the total energy of an isolated Fe dot on 
NaCl(lOO), and the energy of the strain-induced dipolar 
interaction between the dots. 

The experiments were performed in an ultrahigh vac- 
uum (UHV) system with base pressure of lxlO 10 torr. 
The system is equipped with electron beam sources, 
laser molecular beam epitaxy, and an in-situ Omi- 
cron variable-temperature UHV beam deflection atomic 
force microscope (AFM) / scanning tunneling microscope 
(STM) with cooling and heating facilities covering a tem- 
perature range of (13-1500) K. The noncontact mode 
AFM was used to study the surface morphology in this 
work. The NaCl single crystal substrates were cleaved in 
air, then were immediately loaded into the UHV cham- 
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ber and were annealed to 530K for one hour to remove 
surface contamination prior to the experiments. AFM 
images of the cleaved NaCl(OOl) surface show monatomic 
height steps and large terraces with no detectable adsor- 
bates. The iron was evaporated from an Fe wire (5N 
purity) heated by electron beam bombardment at a rate 
of 0.04 ML/min (1ML is equivalent to the nominal sur- 
face atomic density of bcc Fe(110), 1.7xl0 15 atoms/cm 2 ). 
The iron dots can be formed only within a finite temper- 
ature window, since it is known that the morphology of 
the substrate changes dramatically if the temperature is 
above 720K jL3l, whereas at low temperatures, only ran- 
dom clusters or percolated iron films form Jl4| . The dots 
discussed in this Letter are grown at a substrate temper- 
ature of 530K. 
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FIG. 1: Nanometer-scale iron dots grown on NaCl(OOl) with 
a nominal Fe dose of 1.7 ML. (a) A typical NC-AFM image. 
Scan area 500 nmx500 nm. (b) Close up and line profile (150 
nm long) of the dots. Scan area 200 nmx200 nm. The unit 
in the line profile is nm. (c) and (d) Diameter and height 
distributions of the dots shown in (a) and the corresponding 
Gaussian fits. The center and dispersion for the diameter 
distribution are 13.5 nm and 1.1 nm, respectively; and for the 
height distribution are 3.76 nm and 0.36 nm, respectively. 

Nanomctcr-scale dots can be formed directly on the 
NaCl(OOl) surface without a wetting layer, indicating 
that the growth proceeds in the Volmer-Weber mode. 
As shown in the AFM images in Figs. 1 and 2, deposited 
Fe atoms tend to nucleate and form dots even at sub- 
monolayer coverage. The preferred VW growth mode is 
not too surprising, given the large surface free energy dif- 
ference (2.48 J/m 2 for Fe versus 0.18 J/m 2 for NaCl)JI|] 
and the large lattice mismatch between Fe and NaClp6[. 
The most eye-catching characteristic of the dots is their 
size uniformity. Fig. 1 (a) shows a representative AFM 
image of the dots with a nominal Fe thickness of 1.7ML. 
The very uniform dots are randomly distributed on the 



terraces. Fig. 1 (b) is a close up together with a line pro- 
file of the dots. Estimated from the line scan, the dots 
are 3.5 to 4 nm in height, and around 14 nm in diameter 
]l7t . Fig. 1 (c) and (d) are the lateral size and height 
distributions of the dots shown in Fig.l (a) along with 
the corresponding Gaussian fits. Consistent with the line 
profile, the centers of the Gaussian fits of the height and 
diameter distributions of the dots appear at 3.76 nm and 
13.5 nm, respectively. Each distribution has a very nar- 
row width. The dispersions, Ah = ((h 2 ) - (/i} 2 ) 1/2 , of the 
height and diameter distributions are only 0.36 nm and 
1.1 nm, respectively, which are less than 10% of the aver- 
age height and diameter values. Such narrow dispersions 
are comparable with the narrowest dispersions achieved 
in semiconductor quantum dot growth via the SK mode 
0, and are particularly striking because clusters grown 
via the VW mode normally have much broader distribu- 
tions E|. 
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FIG. 2: Typical NC-AFM images and height and diameter 
distributions for dots grown on NaCl(OOl) at different Fe cov- 
erages 8. The image area for each coverage is 150 nmxl50 
nm. The dispersions for the height distributions for the four 
coverages are all 0.36nm; and the dispersions for the diame- 
ter distributions are 0.9nm, l.lnm, l.lnm, and 1.2nm, for the 
four coverages, respectively. 

Another striking observation lies in the coverage de- 
pendence of the dots shown in Fig. 2. As the Fe coverage 
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increases, the dots grow larger and taller, but an opti- 
mal size of the dots is always selected while the narrow- 
ness of their lateral and vertical size distributions is pre- 
served. Our quantitative analysis of the dispersions does 
not yield any significant additional (coverage-dependent) 
broadening of the dot size distributions, consistent with 
the impression drawn by naked eye from the AFM im- 
ages shown in Fig. 2. In Fig. 3, the optimal height and 
lateral size of the dots are plotted as functions of the 
iron coverage, together with the density of the dots. The 
decrease of the dot density with increasing coverage indi- 
cates that the growth of the dots is clearly in the coars- 
ening regime. This observation on the variation of the 
island density makes the persistence of the narrow size 
distributions even more dramatic, because in the coars- 
ening regime, the absolute size distribution of the islands 
is expected to broaden with increasing coverage || . 

Preliminary surface magneto-optical Kerr effect 
(SMOKE) measurements show that these iron dots are 
superparamagnetic in nature. The details of the mag- 
netic investigation, however, are beyond the scope of the 
present paper, which focuses on the formation of the 
highly uniform dots. In order to understand the un- 
derlying formation mechanism of these nanoclusters, it 
is highly desirable to know both their detailed atomic 
structures and their stability. Unfortunately, the iron 
dots formed on NaCl cannot be imaged by STM with 
higher (atomic) resolution. On the other hand, in-situ 
time (from immediately after formation up to 10 days) 
and temperature (from room temperature up to 550 K) 
dependent AFM studies reveal little change in the mor- 
phology of the Fe dots, indicating that the dot arrays on 
the surface are in an energetically favorable configuration 
that is at least metastable. 

At the present, the precise underlying mechanism for 
the formation of these narrow-sized iron dots is not com- 
pletely clear, but the following phenomenological mean- 
field theory, stressing on strain-mediated dot-dot inter- 
actions, is able to reproduce most of the salient features 
observed in the experiment. First, the fact that the pre- 
ferred size of the Fe dots is hardly changed upon anneal- 
ing clearly indicates that the energy per atom E, equiv- 
alently, the chemical potential, of the Fe dot assembly 
must have a local minimum at the observed optimal size 
for a given coverage. This suggests that the qualitative 
form of the energy E follows the generic behavior shown 
in Fig. 4 for a given coverage: When the volume ap- 
proaches zero, E is close to the energy of one adatom; 
when the dot grows very large, E is about the energy per 
atom in a bulk Fe crystal, which is much smaller than 
the energy of one adatom. Between the two limits there 
exists a local minimum, mainly induced by the elastic re- 
laxation caused by the discontinuity of the intrinsic sur- 
face stress tensor at the dot edges [Q. Secondly, the ob- 
servation of different optimal sizes at different coverages 
shown in Fig. 2 suggests that the local energy minimum 
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FIG. 3: The density (a), the peak height (b), and the peak 
diameter (c) of the iron dots as functions of the Fe cover- 
age obtained from the Gaussian fits, in comparison with the 
mean field theory predictions (straight lines). Only a single 
parameter, a — 0.46, is used to fit all the slopes of the three 
independent plots. 
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FIG. 4: Schematic drawing of the energy per atom for the 
Fe dot assembly as a function of the volume V at a given 
coverage, showing the existence of local energy minimum. 

is a function of the iron coverage, implying that dot-dot 
interactions should contribute substantially to the energy 
E. Based on these considerations, we can express £ as a 
function of both the iron coverage 8 and the dot size V 
within the mean field approximation as 

E(9,V) = E (V) + P(6,V) , (1) 

where Eq is the self-energy of an isolated dot and P is 
the dot-dot interaction energy. 

The interaction energy per atom,P(f9, V), is originated 
from misfit-induced strain within the system. By adopt- 
ing the standard dipole-dipole interaction between two 
dots as P(Vi,V 2 ) oc I/iU 2 /r 3 g, we have 

3/2 

P (0,V)= 9v ^, (2) 

where g is a constant related to the interaction strength 
and the effective number of the nearest neighbors of a 
given dot. We note that, on one hand, the long range 
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dot-dot interaction is important enough in influencing 
the energy per atom of the dot assembly, thereby the 
existence of an optimal dot size; on the other hand, it is 
not strong enough to cause the center of mass of a given 
dot to change, thereby the lacking of the spatial order 
among the dots]lS[]. 

The self-energy E (V) of an isolated dot of size V in- 
cludes all possible contributions to the dot, such as the 
strain energy, interface energy, step energy, and kink en- 
ergy, but it does not depend on the total coverage 8. 
In principle, an explicit, but rather complex, expression 
can be obtained from existing studies of strain-mediated 
dot formation in heteroepitaxy, such as Eq. (4.15) of 
Ref.Q. But surprisingly, our detailed analysis based on 
Eq. (4.15)Q shows that no qualitatively correct fitting- 
could be obtained for the range of iron coverages explored 
in the present experiment. In the following, we choose the 
simple expression 



E (V) = bV a 



(3) 



as an alternative phenomenological approach, and show 
that the corresponding energy per atom is capable of ex- 
plaining many of the salient features observed in the ex- 
periment. 

First, by setting the first derivative of the energy per 
atom E(0, V) with respect to the volume V equal to zero, 
we obtain a coverage-dependent local energy minimum, 
at the volume baV a+1 ^ 2 = |# 3 / 2 . This, in turn, yields 
the relation between the density of the dots 8/V and the 
coverage 8 



ba 



-(q+i/2) q 
= -t 
2 



(4) 



With a = 0.46, our theory gives an excellent fit to the 
experimental data as shown in Fig. 3(a). Next, assuming 
that the equilibrium shape of the dots does not change 
with the volume fl^| , we can also obtain the dot height 
(h) and the diameter (D) as functions of the coverage, 
given by h oc 8° 52 and D oc 8 052 , respectively. These 
two relationships also find very good fits to the experi- 
mental data as shown in Figs. 3(b) and 3(c). The fact 
that the exponent a as a single fitting parameter yields 
very good fittings to the slopes of three independent sets 
of experimental observations hints on the physical valid- 
ity of the choice made in Eq. (3). Whereas the precise 
underlying physics for the success of Eq. (3) remains to 
be explored, we suspect that it could be related to the 
unusually large lattice mismatch in the present system 
Jl6| . Finally, we note that the dispersion of the distribu- 
tion of the island heights (or the diameters) shows little 
change with the coverage, whereas our mean field the- 
ory predicts that the dispersion should increase with the 
coverage, Ah oc 8 0A . This discrepancy is not too surpris- 
ing, because within the mean-field theory, the resulting 
optimal size of the islands can be reliably obtained, but 



the exponent of the scaling function around the optimal 
value is typically overestimated, a trend well recognized 
in critical phenomena |2(i[ ]. 

In summary, we have shown that, under proper growth 
conditions, the self-assembly of iron dots on the insu- 
lating surface of NaCl(OOl) leads to the formation of 
nanometer-scale magnetic iron dots with narrow size dis- 
tributions in the absence of a wetting layer. We have 
also demonstrated that the vertical and lateral sizes of 
the dots can both be changed by the iron coverage while 
the narrow size distributions are preserved, even though 
the clustering is clearly in the strong coarsening regime. 
These striking observations have been interpreted suc- 
cessfully with a phenomenological mean-field theory. 
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